Europeanjournalof ScieilCe 



European Journal of Soil Science, July 2014, 65, 573-583 



doi: 10.1111/ejss.l2143 



Ranking factors affecting emissions of GHG from 
incubated agricultural soils 

S. Garci'a-Marco^ S. R. Ravel la''* D . C h a d w i c k''- , A . Va l lejo'=, A. S. Gregory'' & 
L. M. Cardenas'' 

^Departamento de Quimica y Andlisis Agricola, Escuela Universitaria de Ingeiiiena Tecnica Agricola, Universidad Politecnica de Madrid, 
Ciudad Universitaria, 28040, Madrid, Spain , Rothamsted Research, North Wyke, Okehampton, Devon EX20 2SB, UK , Departamento de 
Quimica y Andlisis Agricola, Escuela Tecnica Superior de Ingenieros Agronomos, Universidad Politecnica de Madrid, Ciudad 
Universitaria, 28040, Madrid, Spain , and '^Rothamsted Research, West Common, Harpenden, Hertfordshire AL5 2JQ, UK 



Summary 

Agriculture significantly contributes to global greenhouse gas (GHG) emissions and there is a need to develop 
effective mitigation strategies. The efficacy of methods to reduce GHG fluxes from agricultural soils can be 
affected by a range of interacting management and environmental factors. Uniquely, we used the Taguchi 
experimental design methodology to rank the relative importance of six factors known to affect the emission of 
GHG from soil: nitrate (NOj") addition, carbon quality (labile and non-labile C), soil temperature, water-filled 
pore space (WFPS) and extent of soil compaction. Grassland soil was incubated in jars where selected factors, 
considered at two or three amounts within the experimental range, were combined in an orthogonal array to 
determine the importance and interactions between factors with a Ljg design, comprising 16 experimental units. 
Within this Lj^ design, 216 combinations of the full factorial experimental design were represented. Headspace 
nitrous oxide (NjO), methane (CH4) and carbon dioxide (CO,) concentrations were measured and used to 
calculate fluxes. Results found for the relative influence of factors (WFPS and N03~ addition were the main 
factors affecting NjO fluxes, whilst glucose, NOj" and soil temperature were the main factors affecting CO2 
and CH4 fluxes) were consistent with those already well documented. Interactions between factors were also 
studied and results showed that factors with little individual influence became more influential in combination. 
The proposed methodology offers new possibilities for GHG researchers to study interactions between influential 
factors and address the optimized sets of conditions to reduce GHG emissions in agro-ecosystems, while reducing 
the number of experimental units required compared with conventional experimental procedures that adjust one 
variable at a time. 



Introduction 

Agriculture directly contributes between 5.1 and 6.1 Pg 
CO2 -equivalents (eq) (10-12%) to global greenhouse gas (GHG) 
emissions (Smith et ah, 2007). These emissions are mainly in the 
form of methane (CH4) (up to 3.3 Pg C02-eqyear~'), nitrous oxide 
(NjO) (up to 2.8 Pg C02-eqyear~') and carbon dioxide (CO2) 
(few comparable estimates exist) (Smith et al., 2007). However, 
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agriculture has the potential to mitigate these emissions by using 
existing agricultural technology, with considerable potential in 
adapting cropland management practices (Smith et al, 2007). 

Understanding the processes responsible for GHG emissions from 
agricultural soils, the major controlling factors and which combi- 
nation of factors minimize emissions is critical in the development 
of effective mitigation. It is possible that management and fertilizer 
practices conducted to reduce one GHG could favour conditions for 
production and emission of another. Nitrous oxide can be produced 
by nitrification, incomplete denitrification, nitrifier denitrification 
and/or fungal denitrification, as well as during reduction of nitrate 
(NOj") to ammonium (NH4+) (Kool et al., 2010). Methane is pro- 
duced by methanogenic microorganisms in anaerobic micro-sites 
(Chan & Parkin, 2001) and consumed by soil methanotrophs, espe- 
cially in aerobic micro-sites (McLain & Martens, 2006). Carbon 
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dioxide, the other main GHG generated by soils, is produced by 
autotrophic respiration (roots and associated mycorrhizae) and het- 
erotrophic respiration (soil macro- and micro-fauna) (Hanson et al., 
2000). The link between the N and C cycles in these pathways 
and the influence of different environmental and soil variables 
in each of these processes introduce a large number of factors, 
which affect fluxes from agro-ecosystems (Robertson & Groffman, 
2007). This complexity is especially important for NjO emissions 
because of its large global warming potential. Thus, Snyder et al. 
(2009) listed the most important factors for emissions as (i) 
soil physical and chemical properties, including organic C con- 
tent, mineral N content, pH, texture, drainage, temperature, mois- 
ture content, O2 status, porosity, microbial abundance and activ- 
ity, and (ii) management-related factors, including N application 
rate per fertilizer type, fertilizer application technique, application 
timing, tillage system, irrigation, incorporation of crop residues 
and type of crop. Many of these factors also affect CH4 and CO, 
emissions. 

The individual effects of the major controlling factors on GHG 
emission processes are well known from laboratory-based and field 
studies. However, a lack of knowledge exists on how the interactions 
between those factors affect emissions and how to best optimize 
conditions to mitigate fluxes. 

Laboratory experimentation offers the possibility of a better 
evaluation of treatment effects because it is easier to control 
environmental and soil variables (Schaufler et al, 2010). However, 
conventional experimental approaches that involve altering one 
factor at a time while keeping all other factors constant are difficult 
because of the large number of experimental units this statistical 
approach generates. They are also time consuming and are unable 
to provide the wider range of combinations of interacting factors of 
interest. 

The experimental design (design of experiments, DOE) used in 
the Taguchi methodology is a different approach. It estimates the 
effects of the main factors and the mutual interaction between 
the chosen factors, and determines the optimum experimental 
conditions using only a portion of the total possible combinations of 
a large number of factors (Taguchi, 1987; Rao et al., 2008). While 
traditional experimental design focuses on the average process 
performance characteristics, this approach concentrates on the 
effect of variation on the process characteristics (Phadke, 1989; 
Ross, 1996) and makes the product and or process performance 
insensitive to variation by proper design of factors. 

In this approach, factors are arranged in an orthogonal array (OA) 
to reduce experimental errors and to enhance the efficiency and 
reproducibility of laboratory experiments (Rao etal., 2008). The 
OA properties are such that between each pair of columns each 
combination of levels appears an equal number of times. There are 
several design arrays denoted by a subscript such as L4, Lg, Lg, Ljj, 
Ljg, Ljg, L,7 and and these arrays indicate the main information 
on the extent of the experiment. For example, Ljg has 16 trials, units, 
experimental conditions or combinations of these. Because of the 
orthogonal layout, the effects of the other factors can be balanced 
to give a relative value representing the effects of a level to be 



compared with the other levels of a given factor (Taguchi, 1987). In 
OA experiments, the number of test runs or experimental conditions 
is minimized while keeping the pair- wise balancing property (Byrne 
& Taguchi, 1987). 

The Taguchi method has been mainly and successfully applied 
to improve the quality of manufactured products. To date it has 
been used to limited extent only to determine the optimum process 
parameters in different areas, including microbial fermentations, 
molecular biology, food processing, waste water treatment and 
bioremediation (Rao etal., 2008), and also to study soil erosion 
(Sadeghi et al, 2012). 

To the best of our knowledge, no application of this methodology 
to the study of GHG emissions has been reported until the present 
time. The Taguchi approach could be useful to determine optimum 
conditions to minimize the fluxes, relating these factors and their 
corresponding levels, and reducing the number of experiments 
required. 

The specific objective of this investigation was to apply, for the 
first time, the Taguchi orthogonal experimental design to rank the 
relative importance of different factors on GHG emission processes 
from agricultural soils, and determine which set of conditions are 
required to minimize GHG fluxes. A soil incubation experiment 
(using grassland soil) with five of the most irnportant controlling 
factors (addition of nitrate, carbon (C) quality, soil temperature, 
soil moisture as percentage of water-filled pore space (WFPS) and 
rate of soil compaction) with an OA layout of was performed 
using two or three levels of each factor. Our aim was to identify 
the relative importance of each factor in the emission of GHG from 
agro-ecosystems in order to establish mitigation strategies based on 
agricultural management practices. 

Materials and methods 

Experimental design 

This incubation study was designed by following the DOE method- 
ology and adopting the Taguchi approach. This approach comprises 
the different phases (with various steps) of planning, conducting, 
analysis and implementation. 

In the Taguchi method, performance is measured by the deviation 
of a characteristic from its target value and a loss function [L(j)] is 
developed for the deviation (Ross, 1996), as represented by: 

L(y) = k{y - mf , (1) 

where k denotes the proportionality constant, m represents the target 
value and y is the experimental value obtained for each trial. In the 
case of looking for the quality characteristic of 'smaller is better', 
m is equal to 0 and the loss function can be written as: 

L(y) = ky^, (2) 

and the expected loss function can be represented by: 

E[L{y)\=kEy\ (3) 
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Table 1 Selected factors and their levels (rates) assigned to different soil 
columns used in the Taguchi experiment 





Level 1 


Level 2 


Level 3 


Factors 


Low 


Medium/high 


High 


Nitrate (KN03)/kgN ha"' 


25 


50 


75 


Glucose gCg~' N 


C:N = 0 


C:N = 5 


C:N=10 


Cellulose g C g"' N 


C:N = 0 


C:N = 5 


C:N=10 


Temperature / °C 


15 


25 




WFPS / % 


<80 


>80 




Soil compaction / kPa 


50 


200 





where E{y^) can be estimated from a sample of n as: 

1=1 " 



glucose and cellulose each had three rates or 'levels' assigned (3^ in 
the OA layout), and the soil temperature, percentage of WFPS and 
soil compaction each had two levels assigned (2' in the OA layout) 
(Table 1 and Table SI). These levels were chosen according to the 
expert opinion of the authors, results from previous experiments and 
published literature on known individual controls (Cardenas et al, 
2003; Robertson & Groffman, 2007; Snyder et al., 2009). 

The number of factors and levels determined the selected exper- 
imental design, which was a Ljj OA, therefore comprising 16 
experimental units. In the next step, an experimental matrix was 
designed (h^^-l6 experimental units, with a layout of 2''x3' 
(Level''''"'" X LeveF'"^""')). The total degrees of freedom for the stud- 
ied experimental design was 15 (the number of experimental units 
minus one). Table 2 shows the combination of levels assigned to 
each factor for the 16 experimental units used in the study. The 
experimental conditions were obtained by combining Table 1 and 
the L|f, OA shown in Table 2. 



Phase I: Planning (selection of factors, levels and OA) 

The first step in phase I of the experimental methodology was to 
identify the important factors (operating parameters) to be tested 
whose variations are known to have a critical effect on the overall 
performance of GHG emissions from soils. Six factors, known to 
influence GHG fluxes significantly (Robertson & Groffman, 2007; 
Snyder et al, 2009), were considered: concentration of N03~, C:N 
ratio adding glucose as labile C, C:N ratio adding cellulose as 
non-labile C, soil temperature, soil moisture (as percentage of 
WFPS) and level of soil compaction (Table 1 and Table SI). Nitrate, 



Phase II: Conducting the designed experiments (soil 
incubations and GHG emission measurement) 

Soil was sampled from the top layer (0-10 cm) of a 1-ha field plot 
of permanent grassland on the Rowden experimental platform at 
Rothamsted Research, North Wyke in Devon, UK (50:46: 10°N and 
3:54:05°W). The soil is a clayey pelostagnogley of the Hallsworth 
series in the Soil Survey of England and Wales system (Clayden & 
Hollis, 1984) (a Dystric Gleysol in the FAO, 2006 classification) 
(Table S2). In order to obtain two 'field' replicates, two batches 
of soil were prepared and kept separate; soil samples from five 



Table 2 Orthogonal array Ljg (2^ X 3^^) of DOE (Taguchi methodology) and average of cumulative fluxes (after 16 days) as NtO, CH4, CO2 and total GHG 
as CO2 equivalents, (mean ± SEM, n = 2) 



Cumulative fluxes 



Levels of factors N2O CH4 CO2 CO2 equivalents 

Experimental Soil 

unit Nitrate Glucose Cellulose Temperature WFPS compaction / pgkg~' dry soil / mgkg~' dry soil 



1 


1 


1 


1 


1 


1 


1 


272 ±288 


-131 ±26 


274 ± 26 


356 ±116 


2 


1 


2 


2 


1 


1 


2 


129 ±136 


84 ±95 


394 ± 17 


435 ± 23 


3 


1 


3 


3 


2 


2 


1 


9410 ±887 


373 ±426 


1000 ±327 


3925 ±610 


4 


1 


1 


1 


2 


2 


2 


9154 ±566 


-98 ± 106 


356 ± 14 


3191 ±159 


5 


2 


1 


1 


1 


2 


1 


5864 ±161 


171 ±169 


328 ±45 


72150± 1 


6 


2 


2 


3 


1 


2 


2 


22 075 ±1228 


53 ±202 


335 ±115 


7180 ±262 


7 


2 


3 


2 


2 


1 


1 


378 ±544 


347 ±1 


738 ±85 


862 ±254 


g 


2 


1 


1 


2 


1 


2 


116±14 


173 ±59 


547 ±10 


586±13 


9 


3 


1 


2 


2 


2 


2 


73 949 ±8316 


181±61 


509 ±18 


23 437 ±2595 


10 


3 


2 


1 


2 


2 


1 


35 902 ± 5469 


181 ±63 


935 ±88 


12 069 ±1605 


11 


3 


3 


1 


1 


1 


2 


337 ±94 


180±42 


845 ±78 


953 ±50 


12 


3 


1 


3 


1 


1 


1 


294 ±170 


170±81 


396 ±42 


490 ±9 


13 


1 


1 


3 


2 


1 


2 


234 ±80 


40 ±130 


512 ±84 


585 ±57 


14 


1 


2 


1 


2 


1 


1 


275 ± 246 


127 ±5 


525 ± 1 1 


613±65 


15 


1 


3 


1 


1 


2 


2 


15 684 ±7002 


184 ±37 


300 ±4 


5166±2165 


16 


1 


1 


2 


1 


2 


1 


1879 ±780 


-91 ±4 


259 ± 100 


840 ±141 


Grand average ( 


mean 


± SEM, n = 


16) 








10 997 ±4902 


142 ±28 


516±60 


3927 ± 1528 
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Table 3 Final WFPS (%) and gravimetric moisture content (%) in soils after 
treatment application (mean ± SD, n = 4) 

WFPS / % 





200 kPa 




50kPa 






Field replicate 


>80% 


<80% 


>80% 


<80 


% 


1 


101 ±2 


75 ±2 


92 ±2 


56 ± 


1 


2 


98 ±3 


78 ±2 


93 ±3 


60 ± 


1 




Moisture content / % 








1 


42.8 ±0.7 


34.7 ±0.3 


46.8 ±0.5 


34.9 


±0.2 


2 


42.0 ±0.6 


35.6 ±0.2 


46.7 ±0.5 


35.9 


±0.2 



(field replicate 1) and four (field replicate 2) different points of a 
'W'-shaped transect were taken and pooled in separate containers. 
Both field replicates were immediately transported to the laboratory, 
sieved (< 6 mm) whilst still in a field-moist condition and then 
stored at 2-4° C until required. 

Approximately 270 g of fresh soil was packed into a cylindrical 
plastic core (6.3 cm diameter, 10 cm height) and compressed uniax- 
ially twice to pressures of 50 or 200 kPa with a 2 kN load cell at a 
rate of 1 00 kPa minute"' on a DNIO test frame (Davenport-Nene, 
Wigston, UK) (Gregory etai, 2009) to obtain two batches of 16 
cores per field replicate with two different levels of soil bulk density 
(0.88 ± 0.01 and 0.73 ± 0.01 g dry soil cm-^ respectively, /; = 32). 
Cores were placed in a tray containing approximately 1-2 cm of 
deionized water for 2 days and then half were allowed to drain for 
1 day to bring the soils to two levels of WFPS (< 80 and > 80%). 

As shown by the OA in Table 2, each core received an application 
of N but not all received C. Potassium nitrate (KNO3), dissolved 
in ionized water, was the N source and it was applied at rates 
equivalent to 25, 50 or 75 kg N ha"' (45, 90 and 134mgNkg-' 
dry soil, respectively). For the treatments receiving C, glucose or 
cellulose was applied at a rate equivalent to 125 or 250kgCha~' 
(224 and 448mgCkg~' dry soil, respectively) when the N rate 
was 25kgNha-\ 250 and 500 kg C ha"' (equivalent to 448 
and 896mgCkg~' dry soil, respectively) when the N rate was 
50 kg N ha"', and 375 and 750 kg C ha"' (equivalent to 672 and 
1344mgCkg~' dry soil, respectively) when the N rate was 
75 kgNha"'. There were also some cores that received N but not 
C. This provided amendments with C:N (g : g) ratios of 0, 5 and 
10. The total quantity of water added with C and N amendments 
was 5 ml per core. The final percentage of WFPS of cores is shown 
in Table 3. 

In order to measure emissions of GHGs, 32 packed cores (16 per 
field replicate) were incubated in 1-litre Kilner jars at either 15°C 
(16 cores) or 25°C (16 cores) in constant temperature cabinets for 
16 days to give a total of 32 experimental units. During incubation, 
the lids of the jars were removed to maintain aerobic conditions 
in the headspace and the jars plus soil cores were weighed every 
day and rewetted with deionized water to replace any evaporative 
losses. Measurements of N,0, CH4 and CO2 emissions were made 
on days 0, 1, 2, 3, 6, 8, 10, 13 and 16. After sealing each jar gas 
samples were taken by syringe at 0, 20 and 40 minutes after the 



lids were replaced: the samples were transferred to pre-evacuated 
vials (22 ml) for analysis. Gas concentrations were analysed by 
gas chromatography using a Perkin Elmer Clarus 580 GC with a 
TurboMatrix 110 headspace auto-sampler fitted with an electron 
capture detector (for NjO analysis) and a flame ionization detector 
connected to a methanizer (for CO2 and CH4 analysis) with two 
identical Elite PLOT Q mega-bore capillary columns (Perkin Elmer, 
Shelton, Connecticut, USA). After the last gas sampling, the soil 
cores were analysed for gravimetric soil moisture content, soil min- 
eral N content (N-NH4''' and N-N03~) and water soluble organic 
carbon (WSOC). 

Thirty-two additional cores (with identical treatrnents to those 
used for measuring GHG emissions) were destructively sampled, 
immediately after amendment addition, and analysed for gravi- 
metric soil moisture content, soil mineral N content (N-NH4"'" and 
N-N03~) and WSOC. Ammonium and N03~ were analysed after 
extraction with 2 m KCl (1:2 soihextractant ratio) by automated 
colorimetric determination (Kamphake etal, 1967; Searle, 1984) 
using a segmented flow analyser (Skalar, Skalar SAN^'-^^, Breda, 
the Netherlands). Water soluble organic carbon was extracted with 
Type I water (ASTM classification at 1:5 soihextractant ratio) and 
then analysed with a total organic carbon analyser (Skalar Formacs 
HT CA14 TOC analyser, Breda, ttie Netherlands) 

Phase III: analysis of the experimental data 

Qualitek 4 software for automatic design and analysis of Taguchi 
experiments (Nutek, Inc., Bloomfield Hills, Michigan, USA) was 
used to study the following: (i) assessment of the individual 
influence of each factor on GHG emissions, (ii) analysis of variance 
(anova), (iii) multiple interactions of the selected factors, (iv) 
determination of the optimum factor combination (resulting in 
minimum emissions) and (v) estimation of performance at the 
optimum condition. The anova was the statistical method used to 
assess the relationship between each factor and emission of GHGs. 
The percentage contribution of each significant factor to the total 
variation was calculated as follows (Ross, 1996): 



55' 

P(%)=-^xWO, 



(5) 



where 55^' is the pure sum of squares due to each factor and 552- is 
the total sum of squares and their formulas are: 



55; = SSf - V, X DOFp, 



N ' 



(6) 



(7) 



DOFp represents the degree of freedom associated with each 
factor and is calculated by subtracting one from the number of 
levels of a factor. 55^ is the sum of squares of a specific factor and 
is the error variance. The general formula for calculating the 
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sum of squares of a specific factor (i.e. factor A) for any number of 
levels is: 




where is the sum of observations under level, T is the sum 
of all the observations, n^ is the number of observations under A- 
level, A' the total number of experimental units performed and kj^ is 
the number of levels of factor A. The equation used to calculate 
is: 

F 

V = — = — , (9) 

' DOF^ DOF^ 

where DOF^ represents the degree of freedom associated with error 
and is calculated by subtracting one and DOFp from the total 
number of experimental units performed. 

Interaction between factors was estimated by calculating a sever- 
ity index (SI) of the different factors under study (Roy, 2001). For 
the interaction between factor A and factor C, the generalized for- 
mula for the SI is: 

ABS\{A,C,-A,C,) - iA.C.-AjC.)] X 100 

5/ (%) = ^^-—^ ^^-^ , (10) 

2 X ABS ylargest — smallestl 

where A1C2 is the average effect for factor A at level 1 and factor 
C at level 2 and is calculated by averaging results that contain 
the effects of both A, and and similarly for AjCj, A2C2 and 
A2C]. The strength of the interaction can be measured in terms 
of a numerical quantity that measures the angle between the two 
lines. The SI is formed such that it is 100% when the lines are 
perpendicular and 0% when the lines are parallel. The interaction 
exists when the lines are non-parallel. 

Simple linear correlation analyses were performed by using 
SPSS 15.0 (Chicago, USA) to study the relationships between 
N2O, CH4 or CO2 fluxes and soil properties (percentage 
of WPFS and NO3-, NH4+ and WSOC concentrations) 1 
and 16 days after application of amendment. Total GHGs as 
CO2 equivalents was calculated using the following equation 
(IPCC, 2006): 

C02(eq) = CO2+N2OX 310+ CH^X 21 (II) 

Results 

Relative influence of individual factors 

The GHG emissions were very much dependent on the selected 
factors and their 'levels' (Table 2). The average effects of the factors 
at the assigned levels for each GHG cumulative flux and total GHG 
are shown in Table 4 and Figures S I -S4 in Supporting Information. 
The difference between averages of values at levels 2 and 1 (L2 
minus LI) of each factor is an estimate of the relative influence of 



Table 4 Main effects of individual factors at assigned levels on cumulative 
GHG emissions after 16 days 

Factors Level 1 Level 2 Level 3 L2-L1 



N2O cumulative fluxes / ^ig N kg dry soil 



Nitrate 


4630 


7108 


27 621 


2479 


Glucose 


11470 


14595 


6452 


3125 


Cellulose 


8451 


19084 


8003 


10633 


Temperature 


5817 


16177 


- 


10361 


WFPS 


254 


21740 


- 


21485 


Soil compaction 


6784 


15210 


- 


8426 


CH4 cumulative fluxes / |jg C kg 


dry soil 






Nitrate 


101 


186 


178 


85 


Glucose 


92 


111 


270 


19 


Cellulose 


127 


153 


159 


26 


Temperature 


105 


177 


- 


72 


WFPS 


139 


142 


- 


3 


Soil compaction 


171 


112 


- 


-59 


CO2 cumulative fluxes / mg C kg 


~' dry soil 






Nitrate 


453 


487 


671 


35 


Glucose 


398 


547 


721 


150 


Cellulose 


514 


475 


561 


-39 


Temperature 


391 


640 




249 


WFPS 


526 


503 




-24 


Soil compaction 


557 


475 




-82 


Total GHG as CO2 


equivalents / 


mg C kg~' 


dry soil 




Nitrate 


1889 


2695 


9237 


806 


Glucose 


3954 


5074 


2727 


1120 


Cellulose 


3136 


6394 


3045 


3258 


Temperature 


2196 


5659 




3462 


WFPS 


610 


7245 




6635 


Soil compaction 


2663 


5192 




2529 



the factors on the variability of the GHG's emission process. The 
larger the difference the stronger is the influence. According to this 
criterion our results showed that all factors at level two resulted in an 
increase in N2O emissions, although WFPS was the most important 
factor controlling N,0 emission at an individual level. Table 4 and 
Figure SI also show that N2O cumulative fluxes first increased and 
then decreased when the C:N ratio was largest, regardless of the C 
source used (glucose and cellulose). 

Nitrate addition, glucose at the largest addition and soil tem- 
perature were the most important individual factors increasing 
CH4 emission, although for this gas, soil compaction resulted in 
a decrease of total fluxes (Table 4 and Figure S2). In the case of 
soil respiration (CO2 cumulative fluxes), fluxes were enhanced by, 
in order, soil temperature, glucose addition and nitrate addition, but 
decreased by cellulose, WFPS and soil compaction when level two 
of the factors was applied to the soil cores. However, those fluxes 
also showed an important increase when the largest dose of glucose, 
nitrate and cellulose was added (Figure S3). 

When considering the total GHGs as CO2 equivalents (total of 
N2O, CO2 and CH4), all factors at level two increased the flux, 
whilst the largest C:N ratio decreased them (Figure S4). In general, 
when fluxes of NjO were small, CO2 fluxes contributed more to the 
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Table 5 Example table for analysis of variance (anova) for cumulative N2O fluxes 



Factors 


DOF 


Sum of squares(SS) 


Variance(V) 


F-ratio 


Pure sum (SS') 


Per cent P / % 


Nitrate 


2 


1490141514.56 


745 070757.28 


4470480055 223 


1490141514.56 


32.9 


Glucose 


2 


136 217 580.06 


68108790.03 


40 8612191791 


136217 580.06 


3.0 


Cellulose 


2 


349282 303.06 


174641 151.53 


1 047 905 725 828 


349 282 303.06 


7.7 


Temperature 


1 


429370201.56 


429 370 201.56 


2 576196080706 


429 370201.56 


9.5 


WFPS 


1 


1846528 326.56 


1 846528 326.56 


11078 965 696109 


1846528 326.56 


40.7 


Soil Compaction 


1 


283930925.06 


283 930925.06 


1703 633 704 778 


283 930925.06 


6.3 


Other/eiTor 


6 


-0.001 


-0.001 






0.002 


Total 


15 


4535 470 850.86 











total CO2 equivalent flux, whilst large fluxes of COj equivalents 
were dominated by the N2O contribution. 

Analysis of variance (ANOVA) 

Analysis of variance was used to analyse the results of the experi- 
ment and to determine the percentage of contribution of each indi- 
vidual factor to the variation in emissions. All the factors considered 
had statistically significant effects at the 95% confidence limit. The 
ANOVA table for N2O (Table 5) is included as an example and the 
percentage contribution of each factor on N,0, CH4 and CO2 total 
fluxes, and the sum of the three gases totals in COj equivalents are 
shown in Figures 1 and 2, respectively. 

Optimum levels of the factors for achieving minimum GHG 
emissions 

From the average effects of the factors (Table 4), data were explored 
to determine the conditions that resulted in the minimum GHG 
emissions (Table 6). The expected results at optimum conditions 
were also estimated (Table 6). For N2O, the smallest emissions were 
as indicated in Table 6 and Figure 1(a); WFPS at < 80% and the 
smallest NOj" addition (25 kg ha~') were the dominant factors and 
levels. 

In the case of minimum CH4 emissions (Table 6 and Figure lb), 
the small glucose addition had the major effect (58%), whilst WFPS 
at < 80% had only a small effect (0.02%) at their individual lev- 
els. For minimum CO2 emissions (Table 6 and Figure Ic), the 
most influential factors were glucose and soil temperature, account- 
ing together for 75% of the overall variance of the experimental 
data. For the minimum CO2 equivalents (Table 6 and Figure 2), 
the most influential factors at their individual levels were the per- 
centage of WFPS and NO3" concentration: this was similar to the 
effects on minimum N2O emissions (explaining 74% of the overall 
variance). 

Interaction of factors 

When the interaction of different factors was calculated for each gas 
(Table 7), it was interesting to note that the less important individual 
factors, such as glucose and soil compaction (at their individual 
levels) for N2O and CO2 equivalent had the largest severity index 



percentage (SI) when combined. Hence, these factors had little 
impact on GHG emissions on their own, but had greater impact 
when acting jointly. The larger the SI value the larger the impact on 
the particular process as combined factors. For NjO, glucose with 
cellulose (both at level 2, C:N = 5) and glucose with soil compaction 
(both at level 1, C:N = 0 and 50kPa, respectively) gave SI values 
of 68.8 and 68.6%, respectively. In the case of CH4 emissions 
(Table 7), temperature (at level 1, 15°C) when interacting with 
soil compaction (at 50kPa) or cellulose (at C:N = 5) resulted in 
the largest SI (57.7 and 53.0%, respectively). For CO2 emissions 
(Table 7), soil compaction played a greater role when interacting 
with other factors. Thus the largest SI was found for soil compaction 
(at 50 kPa) in combination with percentage of WFPS at > 80% 
(67.8%), soil temperature at I5°C (48.7%) and glucose at C:N = 0 
(48.0%). When the presence of interactions arnong factors was 
studied for total GHGs as CO2 equivalents, the interacting factor 
pairs found (in order of SI) were the same as those found for 
cumulative NjO fluxes. 

Soil properties 

Although an identical amount of water was added to each soil core 
for the same level of soil moisture, WFPS immediately after water 
application was around 98 and 92%, respectively, in compacted and 
less compacted soils for the large moisture content and about 78% 
(compacted) and 58% (less compacted) for soil cores at the smaller 
moisture content (Table 8). 

The initial mineral N content of the soil before any treatment 
application was 6.0mgN-NO3~ kg"' dry soil and 1.5 mgN-NH4"'" 
kg"' dry soil. The application of amendments resulted in an increase 
in soil mineral N (Table 8), and for NO3" content 1 day after 
application of amendment (1 DAA) this was dependent on the 
amount of KNO3 added. 

The initial WSOC in soil before any treatment application was 
46mgCkg"' dry soil. The addition of glucose or cellulose did not 
enhance WSOC at 1 DAA, or at 16 DAA except in soil cores at the 
greater WFPS (Table 8). 

The cumulative emission of NjO was linearly correlated with 
WFPS at 1 DAA and 16 DAA (respectively, = 0.590 and 0.596, 
P<0.05) and with WSOC at 16 DAA (/^^ = 0.618, P<0.05). 
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(a) Cumulative NjO fluxes 




8% 9% 



(b) Cumulative CH, fluxes 



Compaction 




Figure 1 Percentage contribution of each individual factor to the cumula- 
tive fluxes of N2O (a), CH4 (b) and CO2 (c). 

Significant correlations among cumulative CO2 and CH4 fluxes and 
soil properties were not found. 

Discussion 

The Taguchi methodology provides an excellent opportunity to 
study important interactions between controlling variables of GHG 
production in and emission from soil. Through this increased under- 
standing of soil processes it may allow researchers to develop and 
test strategies to mitigate GHG emissions in agro-ecosystems. The 



Total GHG (CO2 equivalents) 

Compaction 




Figure 2 Percentage contribution of each individual factor to total cumu- 
lative GHG fluxes expressed as CO2 equivalents. 

advantage of the Taguchi methodology is that a large number of 
factors can be studied in the same laboratory experiment with a 
small number of combinations, providing enough information about 
their mutual interactions. In a conventional statistical, full factorial 
design using six different factors with two and three different rates 
or levels for three factors each, it is necessary to have 2^ x 3' (216) 
different combinations. With the DOE design, the number of com- 
binations required to obtain equivalent information as described by 
the OA (Ljj) was 16 as in this experiment. Exploration of the data 
allows researchers to determine the optimized conditions (within 
the given levels used in the experiment) that generate the minimal 
flux, and hence provide essential information on which to base miti- 
gation practices in an effective way. These recommendations would 
then need to be confirmed with field experiments comparing tradi- 
tional management with the best possible options. The results from 
our experiment were consistent with those found in the published 
literature for each GHG where information was available for a lim- 
ited number of factors and/or levels, in other words for one or two 
factors-levels combinations. 

Nitrous oxide 

The factors controlling N2O fluxes were consistent with effects on 
denitrification as the dominant process (instead of nitrification). 
The addition of N03~ as an exogenous source of N, the addition 
of labile C sources, used by heterotrophic microorganisms, and a 
WFPS > 60% (except for three of the present treatments) favour 
denitrification (Davidson, 1991). Our results indicated that the 
percentage of WFPS was the main factor involved (explaining 41% 
of the variance. Figure la) and must be maintained below 80% 
(Table 6) in order to reduce N2O emission: this is in line with 
other findings (Meixner & Yang, 2006; Schaufler et aJ., 2010). 
Although in rain-fed crops, moisture content cannot be controlled 



© 2014 Rothamsted Research. European Journal of Soil Science published by John Wiley & Sons Ltd on behalf of British Society of Soil Science 

European Journal of Soil Science, 65, 573-583 



580 S. Garcia-Marco et al. 



Table 6 Optimum conditions for minimum GHG emissions and tlie contri- 
bution of different factors at optimized conditions. Tlie contribution value 
is the amount of improvement obtained when factor setting is at the ade- 
quate level for minimum emissions. This value is calculated in relation to 
the grand average shown in Table 1 . The expected result at optimum con- 
ditions is an estimate of performance at optimum conditions. This value is 
calculated by adding the contribution from all factors to the grand average 
shown in Table 1 



Factors 



Level description 



Level 



Contribution 



Cumulative N2O fluxes / |jgNkg dry soil 



Nitrate 

Glucose 

Cellulose 

Temperature 

WFPS 

Soil compaction 



Low (25 kgNha"') 
High(C:N=10) 
High(C:N=10) 
Low (15°C) 
<80% 

Low (50kPa) 



-6367 
-4545 
-2994 
-5180 
-10743 
-4212 



Expected result at optimum conditions 
Cumulative CH4 fluxes / |Jg C kg~' dry soil 



-23045 



Nitrate 

Glucose 

Cellulose 

Temperature 

WFPS 

Soil compaction 



Low (25 kgNha"') 
Low(C:N = 0) 
Low(C:N = 0) 
Low (15°C) 
<80% 

High (200 kPa) 



Expected result at optimum conditions 
Cumulative CO2 fluxes / mg C kg~' dry soil 



Nitrate 

Glucose 

Cellulose 

Temperature 

WFPS 

Soil compaction 



Low (25 kg N ha- 1) 
Low(C:N = 0) 
Medium (C:N = 5) 
Low (15°C) 
>80% 

High (200 kPa) 



Expected result at optimum conditions 
Total GHG as CO2 equivalents / mg C kg"' 



dry 



Nitrate 

Glucose 

Cellulose 

Temperature 

WFPS 

Soil compaction 



Low (25 kgNha-l) 

High(C:N=10) 

High(C:N=10) 

Low(15°C) 

<80% 

Low (50kPa) 



1 
1 
1 
1 
1 
2 

-30.0 

1 
1 
2 
1 
2 
2 

114.7 
soil 
1 
3 
3 
1 
1 
1 



-40.4 
-49.4 
-14.5 
-36.1 
-1.4 
-29.5 



-63.3 
-118.1 

-40.9 
-124.4 

-13.1 

-41.1 



-2038 
-1201 
-882 
-1731 
-3317 
-1264 



Expected result at optimum conditions 



-6507 



by farmers, other practices such as fertilizer input and cultivation 
could be synchronized, selecting the best moisture content to reduce 
N2O emission. More possibilities exist in irrigated crops, in which 
irrigation management could be used to control soil moisture 
content. 

The amount of NO3- applied as fertilizer was also an essential 
factor in controlling N2O emission and explained 33% of the vari- 
ance (Figure la). The increase in N emission was not linearly related 
to the N addition (Table 4, Figure SI). Thus, the difference in fluxes 
between LI and L2 (representing an N addition of 25 kg N ha"' ) was 
2.48mgNkg-', whereas that between L2 and L3 (also an increase 
in the N addition of 25 kg N ha"' ) was 20.5 mg N kg"' (Table 4). A 
non-linear relationship has been observed previously for this soil 



Table 7 Estimated interaction of severity index (SI) between two factors 
for cumulative N2O, CH4 and CO2 fluxes after 16 days and total GHG as 
CO2 equivalents 

Factor levels for 

the two factors 

(respectively) 

desirable 

for the optimum 

condition for 



Interacting factor pairs 



Severity index (SI) 

/% 



minimum 
emissions 



Glucose X cellulose 68.8 

Glucose X soil compaction 68.6 

Cellulose X soil compaction 55.9 

Temperature X cellulose 41.9 
Cumulative CH4 fluxes 

Temperature X soil compaction 57.7 

Nitrate X WFPS 56.9 

Cellulose X WFPS 55.6 

Temperature X cellulose 53.0 
Cumulative CO2 fluxes 

WFPS X soil compaction 67.8 

Temperature X soil compaction 48.7 

Glucose X soil compaction 48.0 

Nitrate X WFPS 46.5 
Total GHG as CO2 equivalents 

Glucose X soil compaction 71.1 

Glucose X cellulose 70.4 

Cellulose X soil compaction 56.0 

Temperature X cellulose 41.5 



L2xL2 
LlxLl 
L2xLl 
LlxL2 

LlxLl 
Li xLl 
L2xL2 
LlxL2 

L2xL2 
LlxLl 
LI xLl 
L2xL2 

LlxLl 
L2xL2 
L2xLl 
LI xL2 



type (Cardenas et al, 2010), and these results are consistent with 
our knowledge of enzyme functioning. At large N03~ concentra- 
tions inhibition of N,0 reductase activity occurs (Cardenas et al, 
2003; Bergstermann et a/., 201 1) because of the competitive effect 
of N03~ and NtO as electron acceptors during denitrification. 

The addition of an exogenous C source, such as glucose (soluble 
and easily degradable) or cellulose (non-soluble and gradually 
degradable), together with NO3", produced a contradictory effect 
on emission. The application of a small amount of C relative to 
N (C:N = 5) enhanced emission when compared with no addition 
(level 1, C:N = 0), but mitigated fluxes at a larger application 
(C:N= 10, Table 4), regardless of the quality of added C (glucose 
or cellulose). The availability of C not only supports the activity 
of denitrifiers per se, but also has the indirect effect of causing 
micro-site anaerobiosis, because of increased respiratory demand 
for O2. Increased availability of labile C will favour complete 
denitrification to N2, as also observed by Sanchez-Martin et al. 
(2008) when N was applied with glucose to soil. 

The combination of glucose x cellulose (at level 2 for both 
factors) gave the largest SI values (68.8% in this case; C:N= 10, 
Table 7). An explanation for this could be that glucose could 
prime the use of cellulose (insoluble) by microorganisms. The 
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Table 8 Soil WFPS (%), NO3- (mgN-NOj" kg"' diy soil), NH4+ (mgN-NH4+ kg"' dry soil) and water soluble organic carbon (WSOC, mgC kg"' dry soil) 
contents in the experimental units 1 day after amendment applications (IDAA) and at the end of the experiment (16 DAA) (mean ± SD, n = 2) 



1 DAAWFPS 
Experimental unit / % 


16 DAAWFPS 


1 DAANO3 - 1 6 DAANO3 - 
/ mgNkg"' dry soil 


1 DAANH4-' 
/ mgNkg"' 


- 16 DAANH4+ 
dry soil 


1 DAAWSOC 16 DAAWSOC 
/ mgCkg~'dry soil 


1 


58.1 


±2.8 


59.5 


±3.2 


36.9 


±5.4 


38.2 


± 1.1 


2.6 H 


tO.9 


8.2 ± 


1.5 


40.2 ±7.3 


40.3 ±9.1 


2 


76.7 


±2.6 


72.6 


±1.5 


43.3 


±2.5 


37.0 


± 1.2 


1.9h 


tO.4 


5.7 ± 


1.6 


34.0 ±0.8 


36.1 ±3.0 


3 


91.9 


±0.5 


90.5 


±1.7 


1.1 


±0.8 


0.6 


±0.8 


8.1 H 


1 1.8 


28.9 ± 


11.6 


57. 0± 13.2 


103 ± 12 


4 


98.3 


±4.5 


99.9 


±2.0 


44.1 


±0.7 


2.9 


±2.0 


8.7 H 


t4.5 


44.6 ± 


10.6 


47.6 ±8.7 


70.0 ±8.0 


5 


92.0 


±0.4 


90.8 


± 1.2 


77.4 


±6.5 


24.3 


±5.1 


8.6 H 


t2.2 


18.0± 


0.5 


41.7±13.8 


57.4 ±7.4 


6 


99.1 


±3.1 


101 ± 


3 


79.7 


±0.4 


1.1 


±0.4 


4.1 H 


t2.7 


21.6± 


5.4 


42.4 ±0.2 


74.9 ± 11.4 


7 


58.4 


±4.7 


59.2 


±1.6 


53.6 


±13.0 


58.8 


±2.2 


1.5h 


tO.2 


16.6 ± 


6.9 


36.6 ±0.6 


36.1 ±2.7 


8 


78.5 


±4.1 


74.7 


±3.2 


86.3 


±8.4 


78.6 


±1.3 


4.4 H 


bl.6 


30.0 ± 


7.1 


38.4±17.5 


37.4 ±9.6 


9 


98.9 


±3.0 


101 ± 


2 


109 ± 


11 


0.8 


±0.6 


8.4 H 


t3.3 


41.4± 


8.4 


31.7±6.1 


86.2 ±4.1 


10 


92.9 


±0.6 


96.1 


±3.8 


90.2 


±11.2 


0.1 


±0.2 


8.5 H 


tO.l 


16.8± 


3.6 


44.4 ±5.8 


136 ±24 


11 


78.1 


±3.3 


76.4 


±0.1 


132 ± 


0 


105 ± 


2 


1.2 H 


tO.l 


6.3 ± 


2.9 


21.9±2.3 


27.8 ±2.8 


12 


58.1 


±3.0 


58.6 


±2.6 


116± 


20 


122 ± 


5 


4.5 H 


tl.6 


13.3± 


2.4 


27.7 ± 10.0 


28.9 ± 11.3 


13 


79.1 


±0.3 


75.1 


±2.8 


43.8 


±4.3 


37.6 


±2.5 


3.7 H 


tO.7 


14.9 ± 


2.3 


41.8±3.0 


37.7 ±2.9 


14 


59.6 


±2.8 


58.3 


±1.2 


38.8 


±5.5 


37.3 


±1.6 


1.8H 


tO.7 


13. 8± 


3.5 


26.9±12.1 


34.3 ± 2.6 


15 


99.6 


±3.4 


102 ± 


5 


39.2 


±1.0 


2.7 


±1.8 


3.3 H 


tl.9 


16.1 ± 


6.8 


44.7 ±10.2 


61.8± 10.5 


16 


92.9 


± 1.7 


94.4 


±1.7 


36.5 


±1.1 


0.8 


±1.1 


7.6 H 


tl.5 


16.6 ± 


0.6 


38.1 ±4.4 


81.3±0.4 



combination of both sources of C probably best met the need of 
C by heterotrophic microorganisms in the experimental period and 
maintained a large demand for electrons, which contributed to NjO 
consumption and N2 production. Therefore, strategies based on 
the application of a labile C source in combination with N, such 
as occurs with organic fertilizers or crop residues mixed with N 
fertilizers, must be taken into account and, as occurred in this soil, 
this could have the potential to reduce NjO emission. 

Soil compaction also affected fluxes, increasing from 
6784)igkg~' (less compaction) to 15 210|igkg~' (more com- 
paction) (see Table 4), although the relative importance of this 
factor in relation to others was not great (6.3%). Level 1 com- 
paction represented a static load exerted, for example, by a sheep 
and level 2 that of heavy machinery (Gregory etal, 2007) or 
cow hooves (Scholefield etal, 1985). Increased N,0 emissions 
in compacted soils have been attributed to reduced gas diffusiv- 
ity and air-filled porosity (Skiba etal, 2002). Extrapolating our 
results to soil tillage, it could be expected that soil tillage reduces 
soil compaction (as with level 1) and helps to reduce emissions 
(Menendez etal, 2008). However, optimum conditions for this 
soil could occur when residues with substantial amounts of labile 
C are incorporated by ploughing, because the interaction between 
glucose and soil compaction (SI = 68.6%) and between cellu- 
lose and soil compaction (SI = 55.9%) resulted in less emission 
(Table 7). 

Temperature also increased NjO fluxes. An increase of 10°C 
enhanced the average cumulative NjO fluxes by 2.8 times (Table 4). 
Under denitrifying conditions, a Qjq (the factor by which the 
denitrification rate differs for a temperature interval of 10°C) value 
close to or greater than three has been estimated frequently (Maag 
& Vinther, 1999). 



Methane and carbon dioxide 

Glucose was the main factor to influence CH4 and COj fluxes, 
increasing emissions in line with the rate of C applied. Glucose is 
a labile source that could be easily used by methanogenic bacteria 
to produce CH4 and by other microorganisms producing CO2. The 
other source of C applied in this soil, cellulose, had a small effect 
on emission of CH4 and COj fluxes, probably because of its slow 
rate of degradation in soil. 

Soil temperature is one of the most important regulators of 
ecosystem respiration (Silvola etal., 1996). In our experiment an 
increase of 10°C enhanced soil respiration and CH4 emission by 
1.6 and 1.7 times, respectively (Table 4, Figures S3 and S2). 
This variation is consistent with Qjq values found for biological 
processes, which range from 1 .3 to 3.3 in reviews of soil respiration 
(Tjoelker efa/.,2001). 

Minimization of CH4 requires soil conditions that promote 
methanotrophic rather than methanogenic activity. As expected, 
this was favoured by low temperature, a smaller percentage of 
WFPS (increasing soil aeration) and no application of C to soil 
(Tables 4, 6). Nitrate suppresses CH4 production in anoxic soil 
(Roy & Conrad, 1999), because methanogenesis is inhibited by 
denitrification intermediates (NOj", NjO and NO). However, for 
methanotrophic bacteria, the addition of large amounts of N03~ 
also inhibits CH4 oxidation (Rigler & Zechmeister-Boltenstem, 
1999), probably because of an increase in osmotic potential (Hutsch 
et al., 1996). This may have occurred in our experiment when the 
rate of added N03~ was increased. Our results confirmed that the 
optimum conditions to mitigate CH4 fluxes from this soil were when 
N03~ was applied at the smallest rate (25kgNha~'). Soil com- 
paction and soil moisture only had small effects on CH4 emission, 
which was contrary to our expectation, because more anaerobic sites 
often promote methanogenesis and reduce oxidation. 
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The application of mineral N generally increases CO2 fluxes 
(Iqbal et al., 2009). In our experiment, a non-linear increase was 
found (Ll-L2 = 35mgCkg-' whereas L1-L3 = 184 mg C kg"', 
Table 4 and Figure S3). This increase after N addition was 
probably produced by a priming effect of fertilizer, which 
often promotes a rapid mineralization of soil organic matter 
(Kuzyakov, 2002). 

Total GHG 

Total GHG, expressed as CO2 equivalent, was very dependent 
on NjO emission, except when WFPS was reduced, when COj 
had a larger influence. Therefore, its reduction must be related 
to mitigation of N,0 emission. Other factors not included in this 
study could also affect emissions, for example the use of NH4''" or 
urea-based fertilizers, other C sources with different availability for 
microorganisms, a wider variation of WFPS and the presence of 
plants. 

From a methodological point of view, the use of repacked soil 
cores for this type of study is recommended, because it reduces 
the inherent variability of soil factors and their effect on GHG 
emissions, thus allowing treatment effects to dominate gaseous 
fluxes. This implies that the Taguchi methodology is most suitable 
for laboratory-based experiments, or where uniform soil properties 
would be expected in the field, such as arable-tilled soils or 
homogenous peat soils, but not those from grazed grassland. 

Conclusions 

The Taguchi methodology offers new possibilities to assess com- 
plex interactions between factors known to control soil processes 
and, in the present case, can be used to help establish strategies to 
mitigate GHG emissions from soil. A large number of factors can 
be evaluated in order to quantify their relative effect on emissions, 
but using a reduced number of combinations in comparison with 
those used in an equivalent statistical factor design. In our exper- 
iment, five of the most commonly defined factors affecting emis- 
sions were selected and results were consistent with those found 
in published studies where, in general, conventional experimental 
procedures had been used that adjusted one variable at a time. 

Although other possible factors and levels could have been 
studied, our results indicated that total GHG emissions from 
agricultural soils, expressed as CO2 equivalents, were dependent 
on NjO emissions, where the percentage of WFPS and addition of 
NOj" were the main influencing factors. 

Further studies using other possible factors and rates not included 
in our experiment should be investigated. It would also be of 
value to validate under field conditions the optimized conditions 
that generate minimal and even maximal fluxes obtained with the 
laboratory-based Taguchi methodology. 

Supporting Information 

The following supporting information is available in the online 
version of this article: 



Figure SI. Plot of factor average effects of N2O cumulative fluxes. 
The dashed line represents the grand average, the average of all 
the observations. A line with a small angle or which is horizontal 
has less influence on the variability of the emission process. The 
red circles represent the optimum conditions based on the quality 
characteristic 'smaller is better'. Error bars represent the standard 
error of the mean. 

Figure S2. Plot of factor average effects of CH4 cumulative fluxes. 
The dashed line represents the grand average, the average of all 
the observations. A line with a small angle or which is horizontal 
has less influence on the variability of the emission process. The 
red circles represent the optimum conditions based on the quality 
characteristic 'smaller is better'. Error bars represent the standard 
error of the mean. 

Figure S3. Plot of factor average effects of CO, cumulative fluxes. 
The dashed line represents the grand average, the average of all 
the observations. A line with a small angle or which is horizontal 
has less influence on the variability of the emission process. The 
red circles represent the optimum conditions based on the quality 
characteristic 'smaller is better'. Error bars represent the standard 
error of the mean. 

Figure S4. Plot of factor average effects of total GHG (as CO2 
equivalents). The dashed line represents the grand average, the 
average of all the observations. A line with a small angle or which 
is horizontal has less influence on the variability of the emission 
process. The red circles represent the optimum conditions based on 
the quality characteristic 'smaller is better' . Error bars represent the 
standard error of the mean. 

TableSl. Amounts of nitrate (kgNha"'), glucose and celuUose 
(kg C ha"') added at each C:N ratio. 

TableS2. Physico-chemical properties of the 10 cm of the topsoil 
layer measured by conventional methods. 
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